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ABSTRACT 

The thermal state of the intergalactic medium (IGM) at z < 6 constrains the nature and timing of 
cosmic reionization events, but its inference from the Lyct forest is degenerate with the 3-D structure 
of the IGM on ^ 100 kpc scales, where, analogous to the classical Jeans argument, the pressure of 
the T 10^ K gas supports it against gravity. We simulate the IGM using smoothed particle hy¬ 
drodynamics, and find that, at 2 ; < 6 , the gas density power spectrum does not exhibit the expected 
filtering scale cutoff, because dense gas in collapsed halos dominates the small-scale power mask¬ 
ing pressure smoothing effects. We introduce a new statistic, the real-space Lya flux, Freai, which 
naturally suppresses dense gas, and is thus robust against the poorly understood physics of galaxy 
formation, revealing pressure smoothing in the diffuse IGM. The Freai power spectrum is accurately 
described by a simple fitting function with cutoff at Ap, allowing us to rigorously quantify the pressure 
smoothing scale for the first time: we find Ap = 79 kpc (comoving) at 2 : = 3 for our fiducial thermal 
model. This statistic has the added advantage that it directly relates to observations of correlated 
Lya forest absorption in close quasar pairs, recently proposed as a method to measure the pressure 
smoothing scale. Our results enable one to quantify the pressure smoothing scale in simulations, and 
ask meaningful questions about its dependence on reionization and thermal history. Accordingly, the 
standard description of the IGM in terms of the amplitude Tq and slope 7 of the temperature-density 
relation T = To(p/p)^“^ should be augmented with a third pressure smoothing scale parameter Ap. 
Subject headings: dark ages, reionization, first stars - galaxies: high-redshift - intergalactic medium 
- large-scale structure of the universe - quasars: absorption lines 


1. INTRODUCTION 

The intergalactic medium (IGM) is a rich repository of 
cosmic history. It contains most of the mass of the uni¬ 
verse and records a variety of baryonic and non-baryonic 
processes that occur as the universe evolves (Meiksin 
2009) . This makes the IGM a valuable cosmological lab¬ 
oratory for testing models of structure formation. The 
most readily observable probe of the IGM is the Lya for¬ 
est, the collection of absorption lines seen at redshifts 
2 : ^ 2-5 in spectra of high-redshift quasars. Hydrody- 
namical simulations of structure formation in a AGDM 
universe show that the Lya forest results from the inter¬ 
action between a background of UV photons created by 
active galactic nuclei (AGN) and star-forming galaxies, 
and fluctuations in gas density sourced by gravitational 
instability (Gen et al. 1994; Hernquist et al. 1996; Zhang 
et al. 1997; McDonald et al. 2000; Meiksin et al. 2001; 
Croft et al. 2002 ; Viel et al. 2004; Lukic et al. 2015). This 
picture successfully explains several statistical properties 
of the Lya forest, such as the H l column density distri¬ 
bution, and the line-of-sight power spectrum and proba¬ 
bility density function (PDF) of the transmitted flux. 

The Lya forest can be used to probe the evolution of 
the thermal state of the IGM at redshifts 2 : ^ 2-5. On 
cosmological time scales, the thermal state of the IGM is 

^ Institute of Astronomy and Kavli Institute of Cosmology, 
University of Cambridge, Madingley Road, Cambridge CB3 
OHA, UK; kulkarni@ast.cam.ac.uk 

^ Max Planck Institute for Astronomy, Konigstuhl 17, D-69II7 
Heidelberg, Germany 

^ Heidelberg Institute for Theoretical Studies, Schloss- 
Wolfsbrunnenweg 35, D-69II8 Heidelberg, Germany 

Zentrum fiir Astronomic der Universitat Heidelberg, ARI, 
Monchhofstr. 12-14, D-69I20 Heidelberg, Germany 


expected to evolve through four stages (e.g., Barkana & 
Loeb 2001; Pritchard & Loeb 2012 ): ( 1 ) Between recom¬ 
bination {z ^ 1100) and 2 : = 147, the cosmic gas tem¬ 
perature is coupled to the cosmic microwave background 
(CMB) temperature via Gompton scattering off a small 
fraction of residual free electrons. As a result the gas 
temperature falls as Tgas cx: (1 -h 2 :) from Tgas ^ 10^ K at 
recombination to Tgas ^ 400 K at 2 : = 147. ( 2 ) At lower 
redshifts (z < 147), the free electron fraction is no longer 
sufficient to couple gas and CMB temperatures, and gas 
cools adiabatically as Tgas cx: {l-\-z)^ down to Tgas ^ 2 K 
at ^ = 10. (3) As galaxies start forming at ^ = 10, UV ra¬ 
diation from young, massive stars and possibly a popula¬ 
tion of faint AGN deposits large amount of energy in the 
IGM on time scales of ^ 500 Myr. In this relatively short 
period of time—the “epoch of hydrogen reionization”— 
most of the H l and He l content of the IGM is ionized to 
H ll and He ll respectively. This also heats the IGM via 
photoionization heating so that by 2 : ^ 6 the IGM tem¬ 
perature is Tgas ^ 10^ K (Bolton et al. 2010, 2012). (4) 
Due to the small neutral hydrogen fraction (xhi ^ 10“^) 
in the IGM at the end of the epoch of hydrogen reion¬ 
ization, UV photons can propagate through the IGM to 
large distances without absorption. As a result, a UV 
background is established. The IGM remains in ioniza¬ 
tion equilibrium with this background and begins to cool 
as the Universe expands. However, AGN activity peaks 
at 2 : ^ 3, resulting in an abundance of hard photons that 
ionize He ll to He Hi and likely reheats the IGM to higher 
temperatures. This “epoch of helium reionization” prob¬ 
ably ends at 2 : ^ 2.7 (Shull et al. 2010; Worseck et al. 
2011; Syphers & Shull 2013; Worseck et al. 2014) after 
which the IGM is again expected to cool down to the 
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present epoch. By studying the Lya forest one can con¬ 
strain the thermal state of the IGM providing insights 
into the Universe’s thermal history. 

It is well known that the balance between cooling due 
to Hubble expansion and heating due to the gravita¬ 
tional collapse and photoionization heating give rise to a 
well-defined temperature-density relationship in the IGM 
(Hui & Gnedin 1997)^, which can be written as 

r = r„(y, (1) 

where Tq is the temperature at the mean density p. Im¬ 
mediately after the reionization of H l (^ < 6 ) or He ll 
{z < 3), To is likely to be around 2 x 10^ K and 7 ^ 1 
(Hui & Gnedin 1997). At lower redshifts, Tq decreases 
adiabatically while 7 is expected to increase and asymp¬ 
totically approach a value of 1.62 (Hui & Gnedin 1997). 
The remarkable fact that the mean temperature-density 
relation of the baryons in the IGM can be well approx¬ 
imated by a simple power law is behind the standard 
observational approach of measuring the thermal state 
of the IGM from the Lya forest: statistical properties of 
the forest are measured and compared to hydrodynam- 
ical simulations to constrain Tq and 7 . Statistics con¬ 
sidered for this purpose are line-of-sight power spectrum 
of transmitted flux (Zaldarriaga et al. 2001; Viel et al. 
2009), PDF of wavelet amplitudes in a wavelet decom¬ 
position of the forest (Theuns & Zaroubi 2000; Theuns 
et al. 2002b; Lidz et al. 2010; Garzilli et al. 2012 ), the av¬ 
erage local curvature of the spectrum (Becker et al. 2011; 
Boera et al. 2014), PDF of transmitted flux values (Mc¬ 
Donald et al. 2001; Kim et al. 2007; Bolton et al. 2008; 
Galura et al. 2012 ; Garzilli et al. 2012 ; Lee et al. 2015), 
and the slope of the 6 -parameter distribution (Haehnelt 
& Steinmetz 1998; Theuns et al. 2000; Ricotti et al. 2000; 
Bryan & Machacek 2000; Schaye et al. 2000; McDonald 
et al. 2001; Theuns et al. 2002a; Rudie et al. 2012; Bolton 
et al. 2014). 

However, these measurements rely on the longitudi¬ 
nal, one-dimensional, structure of the Lya forest, which 
is also influenced by the three-dimensional pressure 
smoothing. At large scales and low densities, baryons 
in the IGM are expected to follow the underlying dark 
matter distribution. However, at small scales and high 
densities, baryons experience pressure forces that prevent 
them from tracing the collisionless dark matter. This 
pressure results in an effective three-dimensional smooth¬ 
ing of the baryon distribution relative to the dark matter, 
at a characteristic scale known as the pressure smoothing 
scale. This pressure smoothing affects absorption lines in 
the Lya forest by reducing power in the distribution of 
H I at scales smaller than the pressure smoothing scale. 
This effect, which reduces the small-scale power in the 
longitudinal transmitted flux power spectrum, is degen¬ 
erate with a similar but one-dimensional reduction due to 
thermal broadening (Peeples et al. 2010a,b; Rorai et al. 
2013; Garzilli et al. 2015). 

This degeneracy has been largely ignored by previous 
measurements of the IGM temperature, all of which use 

^ The temperature-density relation is often also called the effec¬ 
tive equation of state of the IGM. The thermodynamic equation of 
state is that of an ideal gas. 


line-of-sight information from the Lya forest (but see 
Becker et al. (2011) and Puchwein et al. (2014)). As 
a result, published constraints on Tq and 7 are confus¬ 
ing and sometimes contradictory. For example, measure¬ 
ments of 7 do not agree with each other; some measure¬ 
ments even suggest an inverted temperature-density rela¬ 
tion (7 < 1 ), which is difficult to obtain in canonical IGM 
models (Bolton et al. 2008; Viel et al. 2009; Galura et al. 
2012 ; Garzilli et al. 2012 ).^ Gonsequently, an important 
goal of IGM studies is to break this degeneracy and bring 
robustness to measurements of the thermal state of the 
IGM. 


The degeneracy can be broken by a suitable transverse 
measurement that is sensitive only to pressure smooth¬ 
ing. Rorai et al. (2013) showed that this can be achieved 
by measuring coherence of Lya forest absorption in close 
quasar pairs. Spectra of close quasar pairs will become 
increasingly coherent as the transverse separation be¬ 
tween the quasars approaches the pressure smoothing 
scale. Rorai et al. (2013) demonstrated that this effect 
can be exploited by measuring the probability distribu¬ 
tion function of phase angle differences between homol¬ 
ogous line-of-sight Fourier modes in quasar pair spectra. 
This new technique has tremendous promise for direct 
measurements of the pressure smoothing scale, and Ro¬ 
rai et al. (2013) showed that applying it to a realistic 
data set should yield measurements with 5% precision. 

The precision with which quasar pair measurements 
will be able to characterize the pressure smoothing scale 
motivates us to better understand pressure smoothing 
in the IGM. Glassically, a simple characterization of the 
pressure smoothing scale is given by (Binney & Tremaine 
2008) 
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in comoving units, where pmo is the average comoving 
matter density, a is the cosmological scale factor, and Cg 
is the sound speed. Thus the pressure smoothing scale is 
a function of the instantaneous gas temperature. Never¬ 
theless, in an expanding universe with an evolving ther¬ 
mal state, the pressure smoothing scale at a given epoch 
is expected to depend on the entire thermal history, be¬ 
cause fluctuations at earlier times expand or fail to col¬ 
lapse depending on the IGM temperature at that epoch. 
In the limit of linear density perturbations, an analyti¬ 
cal estimate of the dependence of the pressure smoothing 
scale on thermal history can be obtained (Gnedin & Hui 
1998). The resultant pressure smoothing scale, termed 
filtering scale, can be written as 
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in comoving units, where pmo is the average comoving 
matter density, 6 is the growth factor of linear pertur¬ 
bations, and Af,J eans is given by Equation (2). Note 
that the filtering scale Af,gh98 is larger than Ap,jeans 
before reionization, and smaller than Ap, jeans (typically 


® This can however be achieved in more exotic models such as 
TeV blazar heating (Chang et al. 2012; Puchwein et al. 2012). 
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Figure 1. Projected density distributions of gas (left) and dark matter (right) at 2 ; = 3 in our fiducial simulation, showing pressure 
smoothing of gas relative to dark matter. The density at each point is an average for a column approximately 5 Mpc/h long. 


by about a factor of 2 ) after reionization. At redshifts 
probed by the Lya forest, typical values of Af,gh 98 5 as¬ 
suming simple models for the IGM thermal state, are 
about a hundred comoving kiloparsecs. With this filter¬ 
ing scale in hand, it is usually assumed that baryonic 
density is given by the dark matter density smoothed 
by a quadratic (Bi & Davidsen 1997) or a Gaussian ker¬ 
nel (Gnedin et al. 2003; Rorai et al. 2013) with width 

Af,GH98- 

However, this description of the pressure smoothing in 
the IGM suffers from two problems. First, scales com¬ 
parable to the pressure smoothing scale, i.e., scales of 
order hundreds of comoving kpc, are highly nonlinear at 
redshifts 2 : ^ 2-5 probed by the Lya forest. Therefore, 
the above linear perturbation theory description of Af 
breaks down at these redshifts. This means that the de¬ 
pendence of the IGM pressure smoothing scale on the 
thermal history is likely more complicated than Equa¬ 
tion (3). Secondly, it is unclear that pressure smooth¬ 
ing of the baryons can be well described as Gaussian 
or quadratic smoothing of underlying dark matter den¬ 
sity, because at scales of ^ 100 comoving kiloparsecs, 
galaxy formation physics could significantly influence the 
baryon distribution. Both of these limitations will affect 
the interpretation of the pressure smoothing scale from 
observations of the Lya forest. The goal of this paper 
is to shed light on the pressure smoothing scale of the 
IGM using high-resolution hydro dynamical simulations. 
We ask if the pressure smoothing in the IGM is well de¬ 
scribed by the approximations in Equations ( 2 ) and (3), 
and present an improved characterization of the pressure 
smoothing scale. 

The structure of this paper is as follows: In Section 2 
we review the quadratic and Gaussian descriptions of 
the IGM and the motivations behind them. We present 
our simulations in Section 3. We analyze our numerical 
results and present a new characterization of the IGM 
pressure smoothing scale in Section 4. Section 5 contains 
a summary of our main conclusions. We work with a 


AGDM cosmological model with = 0.045, Qrn = 0.26, 
Qa = 0.73, h = 0.71, n = 0.96, erg = 0.80 (Hinshaw et al. 
2013). All distances are in comoving units. 

2. ANALYTICAL ESTIMATES 

In this section, we review the linear theory results that 
are traditionally used to describe the pressure smoothing 
scale of the IGM. The IGM is conventionally described 
as an ideal gas with a polytropic equation of state, p = 
Then, in comoving coordinates, linear evolution 
of density perturbations is given by 

^ + 2H^ = 47rC?(pb4 + pM - 

^ + 2H^ = 4nG{p^S^ + pM, (4) 

where S\y{t,k) and are Eourier amplitudes 

of baryon and dark matter density contrasts {6 = 
p/p — 1) respectively and Pb:Pm oc a~^. The ratio 
k)/Sm{t, k) will reveal the pressure smoothing scale 
kp at which pressure support equals gravitational force. 

2.1. Linear evolution for an adiabatie thermal history 

We can ignore pb relative to pm in the gravitational 
source term on the right hand sides of Equations (4). 
The dark matter evolution can then be solved as usual to 
obtain the growing mode b{z) (Peebles 1980). Eor bary¬ 
onic perturbations, a general solution can be obtained at 
2 : > 2 where 6 ^ a, if we assume that the temperature 
at the mean density, Tq evolves as 1/a (Bi et al. 1992). 
This solution can be written as 

S^ipk) ~ l + A^A:^’ ^ ^ 

where we have only retained the growing mode. The 
parameter Aj is defined as 

A^ = 

^ 47rGpm<2^ ’ 


(6) 
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Figure 2. Comparison of the phase-space distributions of gas particles in two of the simulations presented in this paper. (Left panel 
describes the fiducial simulation.) Color scale describes density of points on the plot; most particles are on the temperature-density relation 
shown by the dashed black line. 


which is just the classical Jeans scale of Equation (2). 
Due to Equation (6), for a polytropic IGM, the assump¬ 
tion that To oc 1/a results in a constant Aj at all times. 
Equation (5) describes the behavior of baryons relative 
to the dark matter to linear order in Jb, at redshifts 
2: > 2, under the assumption that Tq oc 1/a. It shows 
that, in linear theory, baryonic density fluctuations are 
damped relative to the dark matter fluctuations by a 
quadratic filter 1 + Aj/c^. At large scales {k Aj), 
baryons follow the dark matter (Jb = ^m) bat at small 
scales {k ^ Aj), baryons are smoothed out by pressure 
support (Jb = ^^m/Aj/c^).^ 


2.2. Linear evolution for an arbitrary thermal history 

However, as discussed in Section 1 above, in general 
the thermal history of the IGM is expected to be dif¬ 
ferent from T oc 1/a. Eor example, T is expected to 
increase during cosmic reionization events, and addition¬ 
ally photoionization heating due to the UV background 
will generally result in different temperature evolution. 
In general, for an arbitrary thermal history T(a) a gen¬ 
eralization of Equation (5) cannot be obtained analyti¬ 
cally. But an approximate result can be derived to first 
order in kf ^ i.e., at large scales (Gnedin & Hui 1998). 
This solution can be written as 


Sb{t,k) 

{^5 


= 1 - 


F5 


(7) 


^ A somewhat more general solution can be obtained for Equa¬ 
tions (4) by assuming To oc which corresponds to a power law 

evolution of Aj (Bi et al. 1992). (The above analysis corresponds 
to a = 1). However, this solution has the same asymptotic behav¬ 
ior as Equation (5). Thus the quadratic filter in Equation (5) is 
a general feature of power law adiabatic thermal histories in the 
linear approximation. 


where 

2 _ 47rGpnio /■*,., KOMn k 
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in which the thermal history enters via Aj and pmo is the 
average matter density at ^ = 0 (recall that Aj is defined 
in Equation (6)). Note that this describes the linear the¬ 
ory evolution of the IGM for an arbitrary thermal history 
only at large scales (/c^Ap ^ 1). A closed-form solution 
of Equations (4) to all orders of kfX^ is not known. How¬ 
ever, by solving Equations (4) numerically, Gnedin et al. 
(2003) show that the full solution is described to a good 
accuracy by a Gaussian of the form 


Jm ifi ^) 


exp (-A:^A|) , 


(9) 


where Ap is defined in Equation (8). Note that this solu¬ 
tion describes the density fields only in the limit of linear 
perturbations. Typically, Equation (9) predicts a pres¬ 
sure smoothing scale of Ap ^ 100 kpc (Gnedin & Hui 
1998). 


2.3. Nonlinear evolution 

Equations (5) and (9) describe pressure smoothing of 
the baryons in the linear perturbation theory regime, 
i.e., in the limit of small density {S < 1). The gener¬ 
alization of these results to the nonlinear limit is un¬ 
known. By arguing that the Lya forest is only sensi¬ 
tive to moderate overdensities (A < 10), Equations (5) 
and (9) have been used in semi-analytical descriptions 
of the forest (Reisenegger & Miralda-Escude 1995; Bi & 
Davidsen 1997; Ghoudhury et al. 2001a, b). Several au¬ 
thors have also used these results to simulate the IGM 
using so-called pseudo-hydrodynamical methods (Petit- 
jean et al. 1995; Gnedin & Hui 1998; Croft et al. 1998). In 
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terms of their predictions for simple cumulative statistics 
of absorption lines in the Lya forest, these methods give 
results that are close to within 10% of hydrodynamical 
simulations (Meiksin & White 2001). 

Nonetheless, the application of Equations (5) and (9) 
to characterize pressure smoothing at the scales and red- 
shifts probed by the Lya forest is problematic. Equa¬ 
tions (5) and (9) typically predict pressure smoothing 
scales of order 100 kpc. At redshifts 2 ; = 2-5, such small 
scales are highly nonlinear, as the variance in the den¬ 
sity distribution at these scales is dominated by collapsed 
structures. Baryons at these densities are affected by 
the physics of galaxy formation, as well as gas outflows 
and inflows. Therefore, in order to understand pressure 
smoothing in the IGM, it would be necessary to some¬ 
how distinguish, at scales of ^ 100 kpc, the low density 
IGM that is likely unaffected by galaxy formation from 
the high density IGM that is affected by it. But where 
to draw the line between these two density regimes is 
not clear. Also, unlike in the linear regime (A < I), 
the pressure smoothing scale Ap is expected to be de¬ 
pendent on the local density in the quasi-linear regime 
(A < 10) probed by the Lya forest, but this dependence 
is unknown, as the classical result of Equation (6) is not 
guaranteed to hold at densities of order 1-10. There¬ 
fore, even if the low density IGM is somehow isolated, 
we are still left with the task of characterizing its pres¬ 
sure smoothing. 

In this paper we address these two issues, and present a 
method of isolating the low-density IGM in a way which 
enables characterization of the pressure smoothing scale 
that not only describes the pressure smoothing in the 
IGM, but also provides an important tool for interpreting 
pressure smoothing scale measurements from the Lya 
forest. 

3. HYDRODYNAMICAL SIMULATIONS 

To understand pressure smoothing in the IGM 
we run high resolution cosmological hydrodynamical 
simulations using the energy- and entropy-conserving 
TreePM smoothed particle hydrodynamical (SPH) code 
p-GADGET-3, which is an updated version of the GAD¬ 
GET code (Springel et al. 2001; Springel 2005). In addi¬ 
tion to the cosmological evolution of baryons and dark 
matter, this code implements photoionization and pho¬ 
toheating of baryons by calculating the equilibrium ion¬ 
ization balance of hydrogen and helium in a optically 
thin UV background, which is taken from an updated 
version of the model of Haardt & Madau (1996) pre¬ 
sented by Dave et al. (1999). Radiative cooling is im¬ 
plemented by taking into account cooling via two-body 
processes such as collisional excitation of H l and He ll, 
collisional ionization of H l. He l, and He ll, recombina¬ 
tion, and Bremsstrahlung (Katz et al. 1996). Likewise, 
P-GADGET- 3 also includes inverse Gompton cooling off 
the GMB (Ikeuchi & Ostriker 1986), which can be an 
important source of cooling at high redshifts. Eor tem¬ 
peratures in the range 10^-10^ K, at redshifts probed by 
the Lya forest, recombination cooling is the dominant 
cooling mechanism (we ignore metal enrichment and its 
effect on cooling rates, which is a good approximation 
for the IGM). 

Baryons with overdensities higher than a few hundred 
times the background density are expected to form galax- 
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Figure 3. Gas (red curve) and dark matter (black curve) density 
power spectra at 2 ; = 3 in our fiducial simulation. Also shown 
are the linearly extrapolated matter power spectrum (dashed grey 
curve), and the quadratic (dashed green curve) and Gaussian (solid 
green curve) filters from Equations (5) and (9). The gas density 
power spectrum does not show a cut-off. 

ies. But as these high density regions are not the subject 
of our study, we simplify their treatment. In the simula¬ 
tions presented here, all gas particles with temperature 
less than 10^ K and overdensity of more than a thou¬ 
sand times the mean baryon density are converted to 
collisionless stars and removed from the hydrodynamical 
calculation (Viel et al. 2004). This speeds up the sim¬ 
ulations, while leaving the low-density IGM unaffected. 
We use the QUICK_LYALPHA flag in P-GADGET-3 for this 
purpose, and in what follows we will also vary the value 
of this threshold to understand its impact on our results. 

Simulations presented in this paper were all done in 
a cubic box of length 10 comoving h“^Mpc on the side. 
Periodic boundary conditions were imposed. The num¬ 
ber of dark matter and gas elements are both initially 
512^. This corresponds to a dark matter particle mass 
of Mdm = 4.5 X 10^ h“^M0 and gas particle mass of 
Mgas = 9.3 X 10^ h ^Mq. This gives us a mean inter¬ 
particle separation of better than 20 comoving h“^kpc, 
which sets our spatial resolution in the lowest density 
regions of the IGM. This resolution is sufficient for simu¬ 
lating the Lya forest (Bolton & Becker 2009; Lukic et al. 
2015). Specifically, at z = 2 it has been found that a mass 
resolution of Mgas = 1.6 x 10 ^ ^“^Mq is sufficient for con¬ 
vergence in the mean Lya fiux and and fiux power spec¬ 
trum in GADGET (Boiton & Becker 2009). At higher red¬ 
shifts, the mass resoiution requirement grows stronger: 
at 2 : = 5 a mass resoiution of about Mgas = 2 x 10^ 
h~^MQ is required. Our mass resoiution meets both 
these requirements (see aiso Bryan et ai. 1999). How¬ 
ever, in their study Boiton & Becker (2009) aiso found 
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Figure 4. Gas density power spectrum at 2 : = 19. The solid red 
curve shows the power spectrum in our fiducial run; solid, dashed, 
and dot-dashed beige curves show power spectra in simulations 
with high To at this redshift. The dotted curve is a Gaussian fit to 
one of the power spectra. Also shown are the linear and nonlinear 
dark matter power spectra. In contrast to Figure 3, the gas power 
spectrum shows a Gaussian cut-off in accordance with Equation (9) 
at this high redshift when most scales in the box are linear. Note 
that in the fiducial simulation the pressure smoothing scale is too 
low to be resolved due to the low temperature at this high redshift. 

that a box size of 10 comoving h~^Mpc is too small for 
convergence in these quantities, giving errors of about 
10% for the flux power spectrum (see also Lukic et al. 
2015). They find that a box size of at least 20 comoving 
h“^Mpc was required ai z = 2. At2; = 5,a box of size 
40 comoving h“^Mpc was required. When the box size is 
small, the error due to the missing large-scale modes and 
the coupling of density perturbation modes in the box to 
modes larger than the box size is large. This results in 
the loss of the overall amplitude of the power spectrum 
of density perturbations (see also Meiksin & White 2004; 
McDonald et al. 2005; Viel & Haehnelt 2006; Tytler et al. 
2009; Lidz et al. 2010; Lukic et al. 2015). However, as we 
will see below, this overall loss of power is not a concern 
for the purpose of characterizing the pressure smoothing 
scale. 

Figure 1 shows distributions of gas and dark matter 
density from the fiducial simulation at 2: = 3 in slices 
approximately 20 comoving kpc/h thick. Gas and dark 
matter distributions are similar on large scales. But on 
small scales, the gas distribution is more diffuse. This dif¬ 
ference is more prominent in quasi-linear structure such 
as filaments and is due to the pressure smoothing that 
we aim to characterize in this paper. 

Apart from the fiducial simulation, we also present re¬ 
sults from four other simulations, which differ from the 
fiducial simulation only in their thermal evolution. In 
these simulations, the photoheating rate of the H l. He l. 
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Figure 5. Evolution of the gas density power spectrum (red 
curves) from redshift 19 to 3 in a simulation in which the photo¬ 
heating rate is enhanced so that the gas temperature is To ~ 10"^ K 
at 2 : = 19. Also shown are the nonlinear (solid black curves) and 
linear (dashed grey curves) dark matter density power spectra at 
each redshift. The gas power spectrum shows a Gaussian cut-off 
as expected from Equation (9) at 2 ; = 19, but this cut-off vanishes 
at low redshifts. 

and He ll in the IGM is enhanced artificially by constant 
factors of 2, 5, 10, and 50 (Becker et al. 2011) resulting in 
different thermal histories. The temperature of the IGM 
at mean density, Tq at 2; = 3 for these simulations is, 

1.4 X 10^ K, 2.4 X 10^ K, 3.6 x 10^ K, and 8.6 x 10^ K re¬ 
spectively. In the fiducial simulation Tq at this redshift is 

9.5 X 10^ K. Figure 2 shows the phase space distribution 
of gas particles at 2; = 3 in the fiducial simulation and in 
the high-temperature simulation with Tq = 8.6 x 10^ K. 
(Photoheating rate in this simulation is enhanced by a 
factor of 50 relative to the fiducial rate.) Both distribu¬ 
tions have familiar features: an overwhelming fraction of 
low density (A < 10) gas particles reside on the power 
law temperature-density relationship of the form given 
by Equation (1) shown by the dashed line. A fraction 
of the highest density gas (A > 100) has cooled in col¬ 
lapsed haloes, and the rest is shock heated to virial tem¬ 
peratures exceeding 10^ K. At this high redshift {z = 3) 
the low-density high-temperature warm hot intergalac- 
tic medium (WHIM; Gen & Ostriker 1999, 2006; Dave 
et al. 2001) is still sparsely populated. The phase space 
diagrams in Figure 2 show a cut-off at A = 10^ for tem¬ 
peratures below 10^ K, which is due to the simple star 
formation recipe described above, which converts these 
cool high-density gas particles to stars. 

For all simulations, we take snapshots of particle po¬ 
sitions, velocity, temperature, and other properties. To 
calculate power spectra, we grid the relevant particles to 
create a density field, using a cloud-in-cell (GIG) scheme, 
taking into account the smoothing lengths of SPH par- 











Pressure smoothing in IGM 


7 


A = 1/fc [cMpc] 



Figure 6. Gas density power spectrum with different star for¬ 
mation thresholds. The solid, dashed, dot-dashed, and dotted red 
curves show the gas density power spectra when the star formation 
threshold is set to 10"^, 10^, 10^ and 10, respectively. The black 
curves show corresponding dark matter power spectra. Increas¬ 
ing the star formation threshold removes high density gas, which 
reduces small scale power. 

tides. After calculating the power spectrum, we decon¬ 
volve the CIC kernel, ignoring small errors due to aliasing 
on the smallest scales (Cui et ah 2008). 

Lya forest spectra are created from the particle dis¬ 
tributions of each snapshot, accounting for the redshift- 
space distortions caused by peculiar velocities and ther¬ 
mal broadening. These are computed by accounting for 
the contributions to density, temperature and velocity of 
SPH particles closer to the sightline than their smooth¬ 
ing length, Doppler shifts due to the bulk velocity, and 
the thermal broadening of the Lya absorption line. This 
procedure results in the Lya flux as a function of wave¬ 
length or equivalently time or distance. Following the 
standard approach, we then rescale the UV background 
intensity so that the mean flux of these extracted spectra 
matches the observed mean flux at the respective redshift 
(Peeples et al. 2010a). We use the Lya flux measure¬ 
ments reported by Faucher-Gigude et al. (2008) for this 
purpose. 

4. CHARACTERIZING THE PRESSURE SMOOTHING 
SCALE 

We now turn to the results of our simulations. Our 
first goal is to compare the simulations with the linear 
theory predictions of Equations (5) and (9). 

4.1. Absence of a Cut-Off in the Gas Density Power 
Spectrum 

Figure 3 shows the dimensionless power spectra of gas 
and dark matter density contrasts ((5 = p/p—1; for gas we 
use p = D^pcr) at z = 3 from our fiducial simulation. The 


dark matter density power spectrum increases steeply 
towards small scales. The amplitude of the gas density 
power spectrum is suppressed relative to the dark matter 
at all scales probed by the simulation. Note that the 
enhancement in the dark matter power relative to the 
linear theory prediction indicates that scales of several 
Mpc are already nonlinear at this redshift. Therefore, 
it is clear that the pressure smoothing scale, which is 
expected to be of order 100 kpc, is already evolving non- 
linearly. 

However, it is striking that in Figure 3 the gas density 
power does not exhibit a cut-off on any scale. This is very 
different from the linear theory expectations discussed in 
Section 2, which suggested a quadratic or Gaussian cut¬ 
off in the gas power spectrum at the pressure smoothing 
scale according to Equations (5) and (9). This discrep¬ 
ancy is highlighted in Eigure 3, where the predictions of 
Equations (5) and (9) are also shown for an arbitrarily 
chosen pressure smoothing scale. The absence of a cut-off 
in the gas power spectrum shows that the linear theory 
formalism of Section 2 does not describe the physics of 
pressure smoothing on non-linear scales at 2; = 3. We 
observe similar behavior, i.e., an absence of a cut-off in 
the gas density power spectrum, in our simulations over 
the complete redshift range probed by the Lya forest 
{z - 2-5). 

To understand this large difference between the gas 
density power spectrum of Eigure 3 and the linear the¬ 
ory predictions of Equations (5) and (9), it is useful to 
consider the gas and dark matter density power spectra 
at redshifts where density perturbations are still linear. 
We therefore look at the power spectra at ^ = 19 in Eig¬ 
ure 4. At this high redshift, density perturbations are 
linear at most scales resolved in our box as indicated by 
the agreement of the dark matter density power spec¬ 
trum with that predicted by linear theory in Eigure 4. 
Scales ^ 100 kpc comparable to the pressure smoothing 
scale are quasi-linear, and the assumptions behind Equa¬ 
tion (9) are thus still valid. We would therefore expect 
the gas power spectrum to exhibit a Gaussian cut-off. 
However, we face a different problem here at z = 19, 
namely that the gas temperature in the fiducial simu¬ 
lation is To = 6 K as reionization does not occur until 
z = 6. At this low temperature, the filtering scale (Equa¬ 
tion (8)) is expected to be ^ 2 kpc from Equation (2), 
much smaller than the smallest scale resolved in the simu¬ 
lation. As a result, there is no clear evidence for a cut-off 
in the gas density power spectrum. Indeed, in Eigure 4, 
the gas and dark matter power spectra from the fidu¬ 
cial simulation agree at almost all scales, verifying our 
expectation that at low temperatures the tiny filtering 
scale ^ 2 kpc implies that gas follows dark matter for all 
scales resolved by the simulation. 

To remedy this, we consider three other simulations for 
which we artificially enhance the H l photoheating rate 
between 2: ^ 6 and 30. Specifically, at redshifts 2: < 6 the 
photoheating rate is equal to that in the fiducial simula¬ 
tion, but it is set to be much higher erg/s) 

at 2; > 6. These simulations with enhanced heating at 
high redshift can be thought of as “early reionization” 
models. Due to this enhanced heating the gas tempera¬ 
ture at 2: = 19 is To ^ 10^ K in these simulations, which 
increases the pressure smoothing scale of Equation (8) 
to be well above the spatial resolution of our simulation 
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Figure 7. Power spectra of the dark matter density (black), gas 
density (red), and F^eal (blue) fields at 2 ; = 3 from our fiducial sim¬ 
ulation. The dashed grey curve shows the matter power spectrum 
prediction for linear growth of perturbations. 

SO that any cut-off will now be visible. Figure 4 shows 
the gas density power spectra at 2; = 19 from these three 
simulations with various enhanced heating rates and gas 
temperatures. The gas power is seen to drop precipi¬ 
tously between scales of 50 and 100 kpc, and the dotted 
curve shows that this cut-off has a Gaussian form, in ac¬ 
cordance with Equation (9). Moreover, as the gas tem¬ 
perature To increases, the cut-off is seen to move towards 
larger scales as expected for the pressure smoothing scale 
from Equation (8). Thus the classical pressure smooth¬ 
ing scale result appears to be valid at very high redshifts, 
when the ^100 kpc pressure smoothing scale is quasi- 
linear, but seems to break down at z ^ 3 when scales 
^ 100 kpc become highly nonlinear. 

This breakdown of the linear theory pressure smooth¬ 
ing scale picture is clearly indicated in Eigure 5, which 
shows the evolution from 2: = 19 to 3 of gas and dark 
matter density power spectra in one of the simulations 
with enhanced photoheating at high redshifts. Compari¬ 
son with the linear theory matter power spectrum (gray 
dashed curves) indicate that scales ^ 100 kpc are already 
highly nonlinear by 2; = 9. Hence one sees that at 2; = 19, 
when the pressure smoothing scale is quasi-linear, the gas 
power spectrum shows a clear Gaussian cut-off at 50 co¬ 
moving kpc, whereas at lower redshifts when this scale 
is highly nonlinear this cut-off vanishes. In what follows, 
we will explain why the cut-off disappears at low redshifts 
and will present a method to reveal it. 

4.2. The Role of High Density Gas 

As we saw above in Eigure 4, the gas density power 
spectrum does not exhibit any cut-off in power at scales 
^ 100 kpc at relatively low redshifts {z ^ 3) in our fidu¬ 


A = 1/fc [cMpc] 



Figure 8. Teal power spectra with different star formation den¬ 
sity thresholds. The threshold, which acts as a crude model of 
galaxy formation, affects the gas density power spectra dramati¬ 
cally but Teal is insensitive to it. 

cial simulation. At these redshifts and small scales, the 
power is completely dominated by highly over dense and 
non-linear collapsed structures, namely halos. Although 
these structures occupy a tiny volume fraction, they nev¬ 
ertheless have a large impact on the power spectrum. 
Recall that galaxy formation is simplified in our fidu¬ 
cial simulation: all gas particles with overdensity greater 
than 10^ and temperature less than 10^ K are converted 
into stars. This results in a stellar density parameter 
= 0.0075 at 2: = 3, which is already factor of three 
higher than the measured value of = 0.0027 at ^ = 0 
(Eukugita & Peebles 2004). We now rerun this simula¬ 
tion, with the same initial conditions, but with modified 
star-formation overdensity thresholds of 10^, 100, and 10. 
At 2: = 3 we find that thresholds of 10^, 10^, 100, and 
10, result in = 0.0066, 0.0075, 0.0092, and 0.0145, 
respectively. As we reduce the star formation density 
threshold, we form more stars but at the same time we 
remove more high density gas from our simulation and 
the power spectrum analysis. 

Eigure 6 shows the resulting gas density power spec¬ 
tra for these different star-formation density thresholds. 
As the density threshold is progressively decreased, and 
more high density gas is removed from the simulation, we 
see that the power on scales smaller than < 100 kpc is 
increasingly attenuated, and the overall amplitude of the 
power on larger scales is also reduced. This behavior is 
not unexpected, and can be understood in a halo model 
picture (Gooray & Sheth 2002). As dense gas above the 
star-formation threshold is removed, we effectively mask 
dense gas in collapsed haloes in the simulation volume. 
The truncation of the small scale power arises from the 
removal of this gas, which can be thought of as sup- 
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Figure 9. Slices of gas (top left) and dark matter (top right) density, and of the -Geal distribution calculated from gas (bottom left) 
and dark matter (bottom right) distributions at 2 ; = 3 in our fiducial simulation. Density at each point is an average on a cubical cell 
approximately 20 ckpc/h on a side. Each slice has a thickness of one cell length. Note that the color scale in panels showing EVeal is 
inverted. The mapping from density to -Geal reduces high density regions to zero. 


pressing the ‘one-halo’ term of the gas power spectrum. 
However, part of the large scale power of the gas density 
also arises from the clustering of these halos, i.e. the ‘two 
halo term’, hence masking haloes results in a reduction of 
large scale power as well. Therefore on large scales, the 
net effect of the star-formation threshold is to introduce 
a linear bias, such that the power spectra for different 
star-formation thresholds in Figure 6 are all parallel to 
each other at small k values. Figure 6 thus shows us that 
the gas density power spectrum is indeed strongly influ¬ 
enced by nonlinear structure. More importantly, as the 
highest density regions are removed, the gas distribution 
begins to show a hint of a cut-off in power. This is seen 
most clearly in the A < 10 curve in Figure 6. 

In other words, to recover the signature of pressure 


smoothing in the moderately overdense IGM (A < 10), 
we need to remove the influence of gas in high-density 
regions, which would otherwise dominate the small-scale 
power. One way of doing this is by applying an ad hoc 
cut-off in gas density as done in Figure 6. However, the 
choice of this density threshold is arbitrary—it would 
need to be made redshift-dependent, as the amount of 
high-density gas that has collapsed to form stars in¬ 
creases at later times, and is thus dependent on the 
galaxy formation model in the simulation, which we have 
crudely parameterized here with a simple star-formation 
density threshold. However, even in more sophisticated 
simulations which attempt to model detailed baryonic 
process in galaxies such as star-formation and feedback 
using a combination of higher resolution and sub-grid 
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Figure 10. Distribution of various quantities along a randomly-chosen line of sight through our fiducial simulation. From top to bottom, 
panels show the dark matter overdensity, gas overdensity, the distribution of F^eal calculated under the assumption that gas follows dark 
matter at all scales (see Section 4.4), and the distribution of F^eal calculated from the gas distribution in the simulation. The bottom panel 
shows the Lyo; fiux along the line of sight, where the effect of longitudinal thermal line-broadening and redshift-space distortions are clearly 
seen. 


modeling (e.g., Schaye et al. 2010; Vogelsberger et al. 
2014; Crain et al. 2015), our findings indicate that the 
the small scale baryon power spectrum will be influenced 
by the details of galaxy formation and the specific sub¬ 
grid implementation. Therefore the challenge is to devise 
a more robust method of removing the influence of high 
density gas on small scale power. 


4.3. Real-Space Flux Reveals Pressure Smoothing in 
Baryons 

In this section we introduce a better method for reveal¬ 
ing the three-dimensional cut-off in the baryon power 
of the moderate overdensity IGM. Given the three- 
dimensional, statistically isotropic gas density field, we 
calculate the real-space Lya flux at each point in the 
box. This quantity, which we denote by Freak is defined 
by 

Freal = exp(-rreal), (10) 


where the real-space Lya absorption optical depth is 
given by 


Teal 


3A3A, 

S7rH{z) 


nm. 


(11) 


Here Aq, = 1216 A is the rest-frame Lya wavelength, 
Aa is the Einstein A coefficient (also written as Aio and 
can be written in terms of an oscillator strength), H{z) 
is the Hubble constant, and nui is the number density 
of neutral hydrogen. The definition in Equation (11) is 
identical to the Gunn-Peterson formula used to compute 
the Lya forest optical depth, except that the convolution 
integral which accounts for the redshift-space effects of 
the peculiar velocity field and thermal line broadening 
has not been included. Indeed, we have chosen to define 
a statistic in real space precisely because these redshift- 
space effects, specifically the line-of-sight smoothing due 
to thermal broadening, are in fact degenerate with the 
three-dimensional real-space Jeans smoothing that we 
aim to study. Hence, F^eai is a three-dimensional, sta¬ 
tistically isotropic field. 

The neutral hydrogen density is given by 


nm 


aR(T)n^ 

Tm 


( 12 ) 


where nu = 0.76pgas is the total hydrogen density, Thi 
is the H I photoionization rate, aR{T) (x is the 

hydrogen recombination rate. Therefore, in a highly ion- 
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ized IGM where the gas lies on a well-defined tempera¬ 
ture density relation T oc (see Figure 2), the neu¬ 

tral hydrogen number density and the optical depth in 
Equation (11) are proportional to a power of the total 
gas density (When all gas is assumed to he 

on the temperature-density relation, Equation (12) con¬ 
stitutes the Eluctuating Gunn-Peterson Approximation 
(EGPA).) Thus, Freai is proportional to the exponential 
of this power of the gas density, and it exponentially 
suppresses the contribution of the high density regions. 
This is of course another way to say that the high density 
regions have large optical depth and therefore result in 
completely saturated hya absorption, i.e., have Freai = 0. 

Freai IS not a directly observable field, but it can never¬ 
theless be constrained by close quasar pair observations 
(Rorai et al. 2013). The phase angle differences between 
homologous line-of-sight Eourier modes in quasar pair 
spectra can be written as 


^12 (^) = COS ^ 


Re 

\FikF2k\ 


(13) 


where Fi(k) and F 2 {k) are Eourier amplitudes of the two 
quasars in a pair. Rorai et al. (2013) showed that the 
probability distribution of these phase differences Ouik) 
is sensitive to a smoothing cut-off in the underlying field. 
Eor modes with wavelength larger than the pair separa¬ 
tion r^, k kj_ ^ l/r_L, the mean value of ^12(^) can 
be written as (Rorai et al. 2013) 


{cos 0{k, r±)) 


Jp dk'k'MrWk'^ - 

jydk'k'Pp^^Jk') 


(14) 

where Jq is the Bessel function of zeroth order and Ppreai 
is the three-dimensional power spectrum of Freai- Thus 
by measuring the phase angle probability distribution 
function in close quasar pairs the Freai can be directly 
constrained. 

Note that in the construction of the Freai held using 
Equation 10, we use the rescaled photoionization rate 
Thi that produces the correct observed mean redshift- 
spaee flux as discussed in Section 3 above. In this ap¬ 
proach, we rescale the UV background intensity so that 
the mean redshift-space flux of these extracted spectra 
matches the observed mean redshift-space flux at the re¬ 
spective redshift (Eaucher-Giguere et al. 2008). When 
this rescaled background is used in calculating the real- 
space flux Freai, hs mean value is not exactly the same as 
the mean redshift-space flux, although the differences are 
very small. Eor instance, in our fiducial simulation, the 
mean redshift-space flux at 2) = 3 is 0.680 after rescaling 
the UV background, while the corresponding mean value 
of Freai, usiug this samc background, is 0.679. Equa¬ 
tions (II) and (12) indicate that rescaling the UV back¬ 
ground to match the observed mean flux values is tan¬ 
tamount to fixing the range of gas densities probed by 
the flux field. As such, the very close agreement between 
the mean of the redshift space F and real-space Freai 
flux indicates that Freai and F probe nearly identical gas 
densities. 

In Eigure 7 we now look at the three-dimensional power 
spectrum of the distribution of Freai at ^ = 3 in our fidu¬ 
cial simulation (recall that unlike the observed Lya flux, 
Freai IS statistically isotropic). Eigure 7 also shows the 


gas and dark matter density power spectra for compar¬ 
ison. The most conspicuous feature of the Freai power 
spectrum in Eigure 7 is that it exhibits a small-scale cut 
off, located at around k = 20 h/Mpc (about A = 35 co¬ 
moving kpc)^. This is reminiscent of the Gaussian filter¬ 
ing scale cut-off that we expect from linear theory. There 
is also an overall reduction in the amplitude of the Freai 
power relative to that of the gas density. This reduced 
power level is similar to that seen when we applied a star 
formation density threshold in Eigure 6, and, which we 
argued can also be understood in the halo model picture. 
Generically, a nonlinear transformation of the density 
field, such as Freai, will have a power spectrum that has 
a different shape than that of the density, but for small k 
where the Eourier amplitudes are small, the transforma¬ 
tion can be linearized such that the power will appear to 
be a linear bias rescaling of that of the original density 
field. 

Thus Freai providcs a natural way of removing the in¬ 
fluence of high density regions in our analysis, and reveals 
the small-scale structure of the IGM—in the form of a 
sharp cut-off in the power at a certain scale—that is pre¬ 
sumably due to pressure smoothing. Eurthermore, this 
definition of Freai is robust against arbitrary choices of 
galaxy formation prescriptions in our simulations. Recall 
that, as discussed in Section 4.2, we can think of varia¬ 
tion in the star formation threshold as variation in our 
(crude) galaxy formation prescriptions. The red curves 
in Eigure 8 illustrate that the gas density power spec¬ 
trum varies dramatically as the star formation density 
thresholds is varied from 10^ to 10. As the threshold is 
decreased, power at small scales is considerably reduced, 
but a clear signature of the pressure smoothing is not ap¬ 
parent from these curves. The removal of clustered high- 
density regions, also reduces the large-scale power as a 
linear-bias rescaling. Although the density power varies 
dramatically depending on the details of how galaxy for¬ 
mation is treated, the blue curves in Eigure 8 show that 
the shape of the Freai power spectrum is invariant to 
these changes. Thus the Freai held as defined in Equa¬ 
tion (10) provides an unambiguous way of characterizing 
the cut-off, valid at all relevant redshifts regardless of the 
galaxy formation details used in the simulation. 

4.4. Visualizing the F^eai Field 

Eigure 9 shows a slice of the Freai field at z = 3 from 
our fiducial simulation in comparison with a slice of the 
gas density field. There is an almost one-to-one cor¬ 
respondence between gas density and Freai- This cor¬ 
respondence results from ignoring redshift-space distor¬ 
tions due to peculiar velocities and thermal broadening. 
The value of Freai is highest (^ 1) in regions with lowest 
gas density, whereas all high density filaments and haloes 
are mapped to Freai = 0. This is exactly the kind of sup¬ 
pression of high densities required to reveal the structure 
of the low-density IGM. Note that unlike in Eigure 1, we 
have chosen a linear color scale in the gas density plot to 
highlight the saturation of Freai values at high densities. 

Eor the sake of illustration, we also apply the transfor- 

® Note that we choose to define the upper-axis as 1/fc instead of 
27T/k because for a Gaussian smoothing the Fourier space cut off 
is inverse of that in the real space. As we see below, the cut off in 
-Geal has a Gaussian form. 
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mation encoded in Freai on the dark matter density field. 
This should be understood as the F^eai held in the limit 
where the baryons directly trace the dark matter density, 
but nevertheless adhere to the nonlinear transformation 
in Equation (12). This is achieved by “gassifying” dark 
matter, i.e., mapping the dark matter density field to a 
pseudo-gas density field through 


f \ 


(15) 


which traces the dark matter at all scales. In other words, 
as the dark matter is cold and pressureless and exhibits 
significant small-scale power, this pseudo-gas distribu¬ 
tion has no pressure smoothing scale. 

The neutral hydrogen number density, uhi, is then cal¬ 
culated by imposing the same temperature-density rela¬ 
tion on these pseudo-gas particles as our fiducial simula¬ 
tion, i.e., by using FGPA. A slice of the resulting 
field calculated in this fashion is shown in the lower right 
panel of Figure 9. Its structure is markedly different from 
that of the F^eai held calculated from the gas distribu¬ 
tion. The F^^^ has a morphology very similar to the dark 
matter density distribution, and because it lacks a pres¬ 
sure smoothing scale cut-off, has a significant amount of 
small-scale structure. In contrast the pressure smooth¬ 
ing of the gas distribution makes the Freai held smoother 
and more diffuse. Although somewhat contrived, F^^^ is 
helpful in better understanding the Freai transformation 
of the density field. 

These points are further illustrated in Figure 10, which 
compares the Freai held along an example line of sight 
through our fiducial simulation to various other quanti¬ 
ties dit z = 3. The upper two panels of this figure show 
the dark matter and gas density field along the line of 
sight. As in Figure 9, we see that the gas density field is 
a smoothed version of the dark matter density field. The 
third and fourth panels of Figure 10 show the F^^^ and 
Freai fields. Calculated from the dark matter and gas den¬ 
sity fields respectively. The mapping to real space flux 
makes the smoothness of the gas density field relative to 
the dark matter density field much more apparent. The 
last panel shows the observed, redshift-space Lya flux, 
F, which includes the effects of peculiar velocities and 
thermal broadening. As a result, the absorption features 
in F are shifted because of peculiar velocities, and their 
widths now reflect the impact of Hubble flow and are 
much smoother because of thermal line broadening. The 
comparison between Freai and F is striking, and illus¬ 
trates that the intrinsic small-scale structure of the gas 
in real space, which is determined by pressure smoothing, 
is completely hidden by redshift space effects. 

Finally, it is instructive to examine the power spec¬ 
trum of the F^^^ field, which is shown as a dashed curve 
in Figure II. As expected, the power spectrum of F^^^ 
shows no evidence for a small-scale cut-off: because the 
dark matter is not smoothed by pressure there is power 
on arbitrarily small scales. The power spectrum at large 
scales is parallel to that for Freai obtained from the gas 
distribution but the amplitude is lower. ^ This compari¬ 
son demonstrates that the cut-off in the power spectrum 


9 The Freai and Fd- fields are non-linear transformations of 
the gas and dark matter density fields respectively. Expressing the 
power spectrum of these non-linear quantities in terms of the origi- 
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Figure 11. Power spectrum of the Feal distribution when cal¬ 
culated under the assumption that gas follows dark matter at all 
scales (dashed curve), compared to the power spectrum of Feal 
calculated from the actual gas distribution in the simulation (solid 
curve). 

of the Freai distribution in our simulations results from 
pressure smoothing. 


4.5. Quantifying the pressure smoothing seale 

We have demonstrated that the power spectrum of the 
Freai field shows a conspicuous small-scale cut-off, which 
we have interpreted as an effect of pressure smoothing. 
In this section we quantify the location of this cut-off, 
and demonstrate that its location behaves as expected 
from the classical Jeans scale argument. 

Inspection of the Freai power spectrum in Figure 7 mo¬ 
tivates a simple fitting function of a power law power 
spectrum cut off by a Gaussian: 

A|(fc) = ^A:"exp ^ (16) 

which has three parameters: A, n, and Figure 12 
shows that Equation (16) is an excellent fit to the power 
spectra of F^eai in out simulations. The fitting function 


nal field requires computing a complicated mode-coupling integral, 
which couples large and small scale power. As such the reduced 
amplitude of relative to F-eal is related to the different shapes 
of the gas and dark matter power spectra. The reduced amplitude 


of is intuitive when one realizes that both Feal and 

span the same domain from zero to unity and have nearly identical 


mean values, whereas the variance of F^^™^ is comparable to that 
of -Feal (although slightly higher because of additional small-scale 
power). The variance of a field is cr^ = l/27r/ A^(fc)dlogA:, which 
is proportional to the area under the curves in Figure 11. Thus 
given the extra small-scale power in Ff^^, in order for F^^^^ and 

-Geal fo have comparable variance, the overall amplitude of Fff^^ 
must be lower. 
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Figure 12. Fits to Freal power spectra from two simulations with 
different temperatures showing that the cut-off in the power spec¬ 
trum is Gaussian. Blue and beige curves show the F^eal power 
spectra from the simulations; black dashed curves are fits. The 
effect of temperature on the cut-off is evident. 

allows US to quantify the pressure smoothing scale in our 
fiducial simulation (Tq = 9.5 x 10^ K) as Aps = l/A^ps = 
79.0 comoving kpc at z = 3. In the simulation in which 
the temperature at the mean density of the IGM, Tq = 
8.6 X 10^ K, is higher than that in the fiducial simulation 
by a factor of ^ 10, is Aps = 191.5 comoving kpc. 

In Figure 13, we compare the F^eai power spectra in 
five simulations with different temperatures ranging from 
T ^ 10^ to ^ 10^. As expected, the dark matter den¬ 
sity power spectra are essentially identical for all tem¬ 
peratures. The gas density power spectrum shows some 
change due to the change in temperature, but these small 
differences would be masked if galaxy formation were 
treated differently (see Figure 6), as we argued in Sec¬ 
tion 4.2. Much more significant changes are seen in the 
Freai powcr spcctrum, where the cut-off scale is observed 
to progressively increase with temperature. Thus Freai 
successfully captures a temperature-dependent smooth¬ 
ing of the gas density field. 

The temperature dependence of the pressure smooth¬ 
ing scale Aps = 1/k^s as determined by fitting Equa¬ 
tion (16) to the Freai powcr spectrum of simulations with 
varying gas temperature is shown in Figure 14. The blue 
and brown curves show the temperature dependence of 
the linear theory Jeans scale, Aj oc \/To, and the lin¬ 
ear theory filtering scale, Ap, defined in Equations (2) 
and (8), respectively. Both of these quantities were de¬ 
rived by measuring the thermal evolution of gas in the 
simulations. As expected at this post-reionization red- 
shift (z = 3), the filtering scale is smaller than the Jeans 
scale, by a factor of ^ 3 in this case (Gnedin & Hui 
1998; Gnedin et al. 2003). We see that the tempera¬ 


Figure 13. Power spectra of the dark matter density, gas density, 
and G-eal fields at 2 ; = 3 in our five simulations with different ther¬ 
mal histories. The dark matter density power spectra are identical 
for all temperatures. The gas density power spectrum shows some 
change due to the change in temperature, but these small differ¬ 
ences would change if galaxy formation were treated differently. 
Much more significant changes are seen in the -Geal power spec¬ 
trum, where the cut-off scale is observed to progressively increase 
with temperature. 

ture dependence of the pressure smoothing scale defined 
by our fitting function in Equation (16) is nearly iden¬ 
tical to the y/T^ dependence of the linear theory Jeans 
and filtering scales, which clearly demonstrates that it 
is probing pressure support in the IGM. Note that the 
slight deviation of our Aps values from the exact \/T de¬ 
pendence is partly due to the fact that simulations with 
different temperatures actually probe somewhat differ¬ 
ent gas densities. This occurs because we always rescale 
the UV background of these simulations to have the same 
observed mean flux value, and by changing the UV back¬ 
ground, one changes the mapping between flux and den¬ 
sity, and hence the range of densities probed by each 
model. This also explains why the pressure smooth¬ 
ing scale is higher than the filtering scale, as the former 
probes pressure smoothing at higher densities than the 
mean density and therefore at higher temperatures than 
To (Schaye 2001). 

This simple fitting function for Treai in Equation (16) 
now allows us to quantify the pressure smoothing of the 
IGM, and make meaningful statements about the pres¬ 
sure smoothing scale for any thermal model and reion¬ 
ization history. 

5. CONCLUSIONS 

In this paper, we presented high-resolution hydrody- 
namical simulations and provided a method to charac¬ 
terize the pressure smoothing scale of the IGM. We pro¬ 
pose a new statistic, which we call the real-space flux. 
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Figure 14. Temperature dependence of the pressure smoothing 
scale at 2 ; = 3 (black curve) in five of our simulations with dif¬ 
ferent values of Tq. The pressure smoothing scale increases with 
increasing Tq. The blue and brown curves show the temperature 
dependence of the linear theory Jeans scale, Aj oc a/Tq? cind the 
linear theory filtering scale, Ap, defined in Equations (2) and (8), 
respectively. The temperature dependence of the pressure smooth¬ 
ing scale is nearly identical to the \/Tb dependence of the linear 
theory Jeans and filtering scales. (The pressure smoothing scale 
is defined by fitting the form in Equation (16) to the finx power 
spectrum, whereas the other two scales in this figure are defined in 
terms of density contrasts. Therefore, for consistency, the pressure 
smoothing scale has been divided by \/2-) 

^reab which is related to observations of correlated Lya 
forest absorption in close quasar pairs and is key to their 
interpretation. Our chief conclusions are: 

• The structure of the IGM in hydro dynamical sim¬ 
ulations is very different from linear theory expec¬ 
tations at redshifts probed by the Lya forest. The 
gas density power spectrum does not exhibit any 
apparent small-scale cut-off, which is the expected 
signature of pressure smoothing in the IGM. 

• This expected cut-off in the gas power spectrum 
is buried under the contribution of dense gas in 
collapsed halos, which dominates the small-scale 
power. As this dense gas is governed by the highly 
uncertain physics of galaxy formation and the cir- 
cumgalactic medium, disentangling the high and 
low density regions of the baryonic field is crucial 
to reveal pressure smoothing. 

• We have introduced a new field, the real-space Lya 
flux, Freai, which naturally suppresses high den¬ 
sity regions associated with galactic halos, and pro¬ 
vides an unambiguous characterization of pressure 
smoothing in the IGM. The power spectrum of Ti-eai 
cuts off at a well-defined scale, which depends on 
gas temperature as expected for the classical Jeans 
scale, making Freai an effective tool for studying 
IGM pressure smoothing. 


• We provide a fitting function that accurately de¬ 
scribes the power spectrum of Freai, and defines the 
pressure smoothing scale as the small-scale cut-off 
of the Freai powcr. This characterization of pres¬ 
sure smoothing enables one to measure the IGM 
pressure smoothing scale from any hydrodynami- 
cal simulation, and thus ask meaningful questions 
about its dependence on the thermal and reioniza¬ 
tion history. Furthermore, it is directly related to 
observational measurements of correlated Lya for¬ 
est absorption in close quasar pair spectra. 

Whereas to date, studies of the thermal state of the 
IGM have focused largely on the amplitude Tq and slope 
7 of the temperature-density relation at a given redshift 
2:, the results of this paper enable one to augment this 
standard description of the IGM with a third parameter 
Aps that describes the pressure smoothing scale, which 
depends on the full thermal and reionization history. 
This characterization of the pressure smoothing scale is 
crucial for interpreting current measurements of the ther¬ 
mal state of the IGM from line-of-sight observations of 
the Lya forest, and is also critical for interpreting mea¬ 
surements of the pressure smoothing scale from quasar 
pair observations. 

Indeed, measurements of the pressure smoothing scale 
using close quasar pairs hold tremendous promise (Ro- 
rai et al. 2013, Rorai et al. in prep.). The sensitivity of 
the pressure smoothing scale to the thermal and reioniza¬ 
tion history of the IGM will be discussed in an upcoming 
paper (Ohorbe et al. in prep.). We will also study the 
requirements for resolving the pressure smoothing scale 
in hydro dynamical simulations, and the convergence cri¬ 
teria for simulating a grid of IGM models with different 
thermal histories (Ohorbe et al. in prep.). Measurements 
of the pressure smoothing scale and comparison to these 
types of theoretical modeling provides a new window for 
understanding the small-scale structure and thermal and 
reionization history of the IGM. The Freai statistic intro¬ 
duced in this paper provides an important new technique 
for relating simulations of the low-density IGM to this 
new observable. 
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